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Abstract 

The SO (4) x U(l) Higgs model on R4 is extended by a F 5 term 
so that the action receives a nonvanishing contribution from the in- 
teractions of 2-instantons and 3-instantons, and can be expressed as 
the inverse of the Laplacian on K4 in terms of the mutual distances of 
the instantons. The one-instanton solutions of both the basic and the 
extended models have been studied in detail numerically. 
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1 Introduction 



Instantons are expected to play a central role in the semiclassical anaysis 
of non-perturbative phenomena in non-linear field theories. Both the Elec- 
troweak and the Strong interactions are described by the dynamics of the 
Yang-Mills(YM) model and hence the SU(2) instantons Q of the relevant 
YM models on IR 4 in each case are expected to be employed. Unfortunately, 
this programme turned out to be less successful than it was hoped. One of 
the main obstacles encountered was the infrared problem arising from large 
instanton effects in the tunneling between topologically distinct vacuua. This 
problem arises directly as a consequence of the scale invariance of the YM 
model on IR 4 , which results in the dependence of the instanton solution on 
an arbitrary scale parameter. The introduction of an absolute scale would 
overcome this problem, and this is the original motivation of the present 
work. 

The most natural way of introducing an absolute scale is by the intro- 
duction of Higgs fields, so that the dimensional parameter representing the 
non-zero vacuum expectation value in the Higgs self- interaction potential 
can set the absolute scale of the would be theory. A related but additional 
feature of the presence of a Higgs field is the exponential localisation of the 
instanton solution, which results in the possibility of constructing a multi- 
instanton field configuaration in which the individual instantons overlap only 
asymptotically. This is what would enable the construction of a dilute gas of 
instantons. Furthermore with the appropriate asymptotic behaviour of these 
instantons, it may become possible to construct a dilute gas of instantons 
which interact non- trivially to yield a Coulomb gas in the appropriate di- 
mensions. By a Coulomb gas we understand a gas where the contribution to 
the action coming from the interactions is given by the inverse of the Lapla- 
cian in terms of the mutual distances of the constituents of the gas. It is 
our aim in the present work to propose a gauged Higgs model on IR 4 which 
acheives the objectives just stated. 

The construction of such a Coulomb gas for a Higgs model on IR 3 was per- 
formed long ago by Polyakovj^]. The Higgs model there [Q] consisted of the 
£77(2) YM field interacting with an su(2) valued Higgs field and the usual 
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symmetry breaking Higgs potential. Firstly, the well known monopolef^] so- 
lution of this model was taken to be the exponentially localised instantons 
in 2 + 1 dimensions. Secondly, a Coulomb gas of these instantons was con- 
structed using the asymptotic fields of the monopole[|J. Our plan in the 
present work is to propose a gauged Higgs model which supports instanton 
solutions with the requisite asymptotic properties capable of describing a 
Coulomb gas on (R 4 . We shall restrict ourselves to this first task here, and 
defer the second technically complicated task of constructing the resulting 
dilute gas action to a future work. Before proceding we make a remark aimed 
at putting the task at hand in perspective: Polyakov's construction [FJ] of the 
Coulomb gas of instantons on IR 3 is the 3 dimensional analogue of the 2 di- 
mensional Coulomb gas of instantons contructed previously by Berezinsky 
and by Kosterlitz and Thouless[[|] employing the 0(2) model on IR 2 , while 
the present work proposes the corresponding instanton field configuration on 
IR 4 . 

Our sole task in this paper is to find the instanton solution of a particular 
Higgs model satisfying the above stated criteria. We shall not attempt to 
compute the action of the corresponding Coulomb gas. Our plan is based 
on the 3 dimensional example 0, namely to construct an appropriate gauged 
Higgs model in contrast with the non-gauged model employed in the 2 dimen- 
sional example^]. Section 2 will be concerned with the instanton solution 
of the basic SO (A) x Z7 (1) Higgs model proposed in Ref.0. We will make 
a detailed numerical study of this solution and find out that inspite of the 
instanton in question being exponentially localised, the asymptotic proper- 
ties do not permit the construction of a dilute gas whose inter-instanton 
interactions support a Coulomb gas. Then in Section3 we shall propose 
two extended versions of the basic model and will verify the existence of 
unit topological charge instantons by numerical techniques. These extended 
models are designed specially to support a non-trivial dilute Coulomb gas 
on R4. Finally in Section 4 we will summarise our results and give a brief 
discussion of the outlook for the extended models. 
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2 The basic model 



First we briefly describe the SO (4) x £7(1) model introduced in Ref. |jj and 
then precede to subsections 2.1, 2.2 and 2.3 where we give the asymptotic 
solutions, the numerical solutions and, an analysis of the inter-instanton 
interactions, respectively. 

The basic model on R4 is described by the antihermitian SO(A) x U(l) 
gauge connection A^ — + 275^ and the antihermitian Higgs multiplet 
$ = 757a0a, A^ being the so(4) connection, a M the U{1) connection and 
(p a a 4-vector real Higgs field. The model, which is derived from the 8- 
dimensional member of the scale invariant Yang- Mills hierarchy on FR 4 x S 4 
by dimensional descent |7]], has the following Lagrange density 

U = TrfF 2 ^ + 4A!{iV, D p ]&} 2 - 18\ 2 ({( V 2 + $ 2 ), F, v } - [D^, D u] $]) 2 
- 54A 3 {(7? 2 + $ 2 ), Zyf} 2 + 54A 4 (r7 2 + $ 2 ) 4 ] (1) 

In ([I]) and everywhere below, the brackets [, ] and {, } denote commutators 
and anticommutators respectively, and square brackets around indices signify 
total antisymmetrisation. Also the curvature F and the covariant derivative 
D pertain to the 5*0(4) x U (1) connection A = A+ij^a, while below we shall 
use F and D pertaining to the SO (4) connection A. The 4-form curvature 
field -F^pcr, in this notaion, is defined as F^^ = {F^ u , F pa ]}. The last term 
in ([!]) multiplying A4 is the Higgs symmetry breaking potential. It is clear 
that any one of the four dimensionless coupling constants A a , a = 1,2,3,4, 
which we shall denote with a vector notation A, may be set equal to 1 by 
rescaling the $ field, but we do not do that here. 

A physically very important property of the system (|l|) is, that at high 
momenta or small distances when the magnitude of the Higgs field is negligi- 
ble with respect to the dimensional constant 77, the system is dominated by 
the term with the highest power of 77. This happens to be 

A** - Tr[~K - HD^) 2 + ±tf + $ 2 ) 2 ], (2) 
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namely the usual Yang-Mills-Maxwell-Higgs system which appears as the 
coefficient of the rf term in (]!]), which is its dominant contribution at high 
energy leading to a perturbative action as expected. Thus the model ([I]) can 
be regarded as the corresponding effective action density 

There are two other physically relevant properties of the system ([I]). The 
first is that the gauge group SO (4) leads to a chirally symmetric dynamics, 
hence does not fit in with Electroweak theory. The other is that at low 
momentum or large distances, the gauge field dynamics is dominated by the 
4-form curvature whose propagation properties again do not fit in with 
Electroweak dynamics. We must therefore conclude that the model ([l|) is not 
suited to be a prototype for an Electroweak theory, but it is consistent with 
the features Strong interaction dynamics both at high and low momenta. 
In this respect, the U(l) connection a M in the gauge field multiplet A = 
{A^ + i^cift) is irrelevant and can consistently be suppressed. We have left 
a p in the definition of £4 in (JJ), because that is the most general || model 
descending from the 4p dimensional member of the Yang-Mills hierarchy that 
can support instantons, and also because when we restrict to the spherically 
symmetric configuration below the imposition of this symmetry will eliminate 
this U(l) field. 

In view of the above reasoning, we will study the properties of this model 
in the context of Strong interaction dynamics at low momentum, and in 
particular will attempt to construct a dilute gas of instantons that satifies 
the properties of a Coulomb gas in R 4 , analogously to Polyakov's work0| in 
R3. As the title states however, this programme will not be completed but 
we will restrict henceforth to the study of the classical instanton solutions of 
this model and its extensions. 

When the values of A a respect the following restrictions: A a > 1 for each a, 
then the action density ([[]) is bounded from below by the topological charge 
density g = d^fl^, expressed as the divergence of the Chern-Simons form il M 
given by|| 

1 2 * * 

^ = --^ £ ^ Tr l5[v A v {F p(T - -A p A a ) 

- \rf<S>{F^ - U({(v 2 + &),F pa } + D a] $)D u $}. (3) 
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The surface integral of this density is the topological charge N, and the 
spherically symmetric solutions with which we will be concerned in this paper 
yield N = 1 for this integral. 

In the limit where X a — 1 for all a, or A = (1, 1, 1, l),this bound is 
saturated and the action is equal to the topological charge N. If however 
any of the A a take positive values smaller than 1, then the action is not 
bounded any more by N but by \ m inN, where \ m in is the smallest of those 
A a which are smaller than 1. This can be shown using similar arguments to 
those used in Ref. 0] in the case of the Abelian Higgs model. When A a < 
for any a, we lose the topological lower bound. 

In the case where the topological lower bound is saturated we have the 
following self-duality, or Bogomol'nyi, equations: 

^W7sV - (V 2 + $ 2 ) 2 = (4) 

' 75 {F [pa , D^} - {(r/ 2 + $ 2 ), = (5) 



-e 



^75({(^ + $ 2 ),^} + A^ D a] &)-({(7 1 2 + $ 2 ),F, U } + D [ ^b u] <5>)=0. 

(6) 

The system (H) shares an important property with the SkyrmefjOJ model 
on IR 3 , namely that all but the second power of the derivative of any field 
component are excluded. This has been a criterion in the construction of 
(lip. Like the Skyrme model also, the self-duality equations saturating the 
topological bound are overdeterminedp| [p"T|j . Accordingly the only solutions 
of the equations (§)-© are the trivial constant solutions we shall see below, 
and any non-trivial solutions we find will satisfy the second order Euler- 
Lagrange equations and not even for Aj = 1 for all i. 

The most important property we shall require of the instanton solution 
is the large r limit of the magnitude of the Higgs field |$| = V — Tr$ 2 

lim |$| = 77 (7) 

in terms of rj, the vacuum expectation value of the Higgs field. Requiring this 
asymptotic condition, we shall find out that the corresponding behaviours of 
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the 50(4) gauge connection and the full Higgs multiplet are 

lim = -g~ x d^ g , lim $ = rjg, (8) 

in terms of the SO (4) group element 75 j^x^, with x = -. We now adapt the 
arguments in Ref.[|l^] pertaining to the SU(2) YM instanton to the solutions 
satisfying (||). In the temporal gauge A Q = each component of A^ and $ 
become time independent, thus enabling the identification of our solutions 
as instantons. 

The solutions we shall seek below are restricted to the spherically sym- 
metric restriction of the Euler-Lagrange equations, namely the unit topo- 
logical charge instantons. This restriction is quite adequate since instanton 
configurations of all chrges can be attained by considering the collection of 
arbitrarily many N = 1 intantons, since as we shall see below these over- 
lap only asymptotically. The reason is the presence of the Higgs field which 
introduces the absolute scale with respect to which the instantons will be 
localised. Indeed, we shall be able to acheive an exponential decay yielding 
this localisation, similar to the instantons of the 4 model [|13|] on IR 2 and 
'compact electrodynamics' [31 |i~3| on IR 3 . 



2.1 Asymptotic solutions 

Since we are unable to solve analytically the Euler-Lagrange equations of 
the one-dimensional subsystem of ([!]) arising from the imposition of spher- 
ical symmetry, we will give explicit solutions in the r >> 1 and r << 1 
asymptotic regions only. 

Under the imposition of spherical symmetry on (R 4 , the Abelian field f^ iu = 
d^a^ defining the system (|I|) vanishes and we are left with only the S0(4) 
field Afj, = A^, for which we employ the following spherically symmetric 
Ansatz 



1 + /(r) 

= l^uXu, $ = h(r) 757m^> ( 9 ) 
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where f(r) is a dimensionless function and h(r) has the same dimensions as 
V- lnv — —\[lpL,lv] is the (Dirac) spinor representation of SO (A). 

Substituting (|9|) into ([[]), performing the angular integrations and with a 
convenient rescaling, we have the following one- dimensional subsystem 



^ = 4(i-m r 2 + ^rv-^ 

+y([(l-f 2 )h]r) 2 + \2r([(v 2 -h 2 )f] r ) 2 

+^(/Ml-/ 2 )) 2 + 4A 3 r(/M^-^)) 2 

+ y ((1 " fW " h 2 ) + 2A 2 ) 2 + ^( V 2 - h 2 )\ (10) 

having suppressed the overall constant factor coming from the angular vol- 
ume. The subscript r in (|10D denotes the differentiation f r = %■ 

Using the notation 5fL = |j — we express the two Euler-Lagrange 

equations of fllQl) . The equation 5fL = arising from the arbitrary variations 
5f is 

S f L 4 = ^(1 - f) (2// f 2 + 2(i _ f)f r - (1 _ f )/ rr ) 

(t(/ 2 - ^ " ^[2/^ + 2// r ^ + 4//A + (/ 2 - l)/i rr ]) 
-2A 2 (/i 2 - ?? 2 ) ([(/i 2 - ?? 2 )/] r + r[2/i 2 / + 2hh rr f + Ahh r f r + (ft 2 - V 2 )f rr }) 

+ 6 ~^fh 2 (f - 1)(3/ 2 - 1) + 8X 3 rfh 2 (h 2 - V 2 ) 2 

+ ^/(3/^ 2 - V 2 ) {{h 2 - v 2 )(f 2 " 1) + V 2 h 2 ) , (11) 

and using the similar notation ShL = |£ — -§p(-§^-), the equation = 
arising from 5/i is given by 

5 h L A = -^r 2 (h 2 - T, 2 ) (2rhh 2 + (h 2 - V 2 )(rh rr + 3h r )) 

7 



+ %/ 2 - 1) ([(/ 2 - l)h) r - r[2f 2 h + 2ff rr h + AffMf ~ ±)hr 



-A\ 2 hf [\{K 2 - V 2 )f} r + r[2ti 2 J + 2hh rr f + Ahh r f r + (K 2 - n 2 )f rr 
+ ^(/ 2 - Iffh + 24X 3 rf 2 h(h 2 - n 2 )(3h 2 - V 2 ) 



+ ^h(3f 2 - 1) ((h 2 - V 2 )(f 2 - 1) + 2h 2 f) + 2A 4 r 3 (/i 2 - V 2 ) 3 h (12) 

Since the system is a non Abelian Higgs model like the models em- 
ployed in Refs.|§, we expect the solution to have asymptotic behaviour sim- 
ilar to that of the monopole||. Accordingly we seek solutions with the fol- 
lowing properties 

lim/(r) = -1 SS>/( r ) = ( 13 ) 

lim/i(r) = lim h(r) = r\. (14) 

The first members of both ( JT3 ) and (FJ) guarantee that the solution is regular 
at r = 0, while the second members satisfy criteria of finite action. Moreover, 
the second members of these equations result in the asymptotic behaviour 
anticipated in fl8|), which is essential for the interpretation of the finite action 
topologically stable solution as the instanton of the dynamical model. 

Expanding around their asymptotic values f(r) = + F(r) and h(r) = 
r] — H{r) in the region r >> 1 in terms of the small functions F(r) and 
H(r), and retaining only linear terms in these functions, ( P]) and (|T2| ) reduce 
respectively to 

p 2 F pp - 3pF p - 3p 2 F = (15) 
o 2 E au - oE a - 8a 2 H = 0, (16) 
in which we have used the following two dimensionless rescalings of r: p = 



2Ai 

The solutions of ( |I~5D and ([U]) respectively are expessed in terms of modified 
Bessel functions and lead to the following exponentially decaying solutions 



/( r ) = A^V^y^r) (17) 
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h(r) =7/(1 



(18) 



in the r >> 1 region. 



In the r << 1 region, we tried a solution in powers of the rescaled radial 
variable p = rjr and found 



/(p)= V + o(p 4 ) (19) 
%) = Sp + (p 3 ) (20) 

where the dimensionless constant A and the constant B with the dimensions 
of r\ are arbitrary and will be determined by the numerical integrations imme- 
diately below. Note that the behaviours ([H]) and (|20|) lead to fields (A^, $) 
(||) that are regular and differentiable at r = 0. 



2.2 Numerical solutions 



We have integrated the equations <5/L 4 = and £>/iL 4 = given respec- 
tively by ([11]) and (|T2| ) numerically starting with the functions ([19]) and (f20|) 
near r = and choosing the numerical values of the constants A and B 
so that f(r) and h(r) reach their asymptotic values given by the second 
members of ([13]) and ([14]). We have performed this integration for the fol- 
lowing values of the dimensionless coupling constants Am = (Ai, A2, A3, A4), 
with A« = (0.5,0.5,0.5,0.5), A (2) = (0.8,0.8,0.8,0.8), A (3) = (1,1,1,1), 
and A 1 - 4 -* = (1.2, 1.2, 1.2, 1.2). The numerical values of the pair of constants 
{A, B} are fixed by each of these the numerical integrations, and are listed 
respectively in the first and second columns of Table 1. The profiles of the 
functions /(r) are given in Figure 1, and those of h(r) are given in Figure 2. 
We do not exhibit the profiles of the action densities corresponding to these 
solutions since they are all ball shaped as expected of spherically symmetric 
lumps. The values of the action integrals S , (A ( -*- ) ), i = 1,2,3,4 for each of 
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these solutions are listed in the third column of Table 1. We defer a detailed 
discussion of these quantitative results to Section 4. 



The above integrations demonstrate the existence of the spherically sym- 
metric solutions of the basic model ||, but we find it interesting to pursue 
our numerical studies somewhat further because of an unusual feature of the 
model at hand. This concerns the overdeterminedjl I ] nature of the Bogo- 
mol'nyi equations (0)-(|§), a feature shared with with the Skyrme model pi? 



on IR 3 . While in the usual Skyrme model the overdetermined self-duality 
equations are parametrised only by one dimensional constant which sets the 
absolute scale, here in addition to the corresponding dimensional constant 
r\ we have the three independent components of the four dimensionless cou- 
pling constants A. It is therefore interesting to investigate numerically the 
quantitative departure from the 'minimal' configuration A = (1,1,1,1). 



We know from Refs.|| |TT| that are only satisfied by the constant 

asymptotic values of /(r) and h(r), given by ([13|) and (|I4"D , and therefore that 
the stress tensor 7^„ corresponding to these solutions does not vanish. It is 
interesting to see how the component 7^ 4 of the stress tensor behaves in detail, 
quantitatively. Now the component T44 of the stress tensor corresponing 
to the spherically symmetric sub-system equivalent to the one- dimensional 
Lagrangian fllPf) is 



T 44 = (^(1 - pff r - - /> 2 ) 4 )+(^rv - h-fhi - ^fh\\ - fr 

+ (±[(l-f)hf r -4\ 1 rfh\ V >-h 2 )) 



+ A 2 r[(^ - tf)fTr - f[0- - f)(v 2 - h 2 ) + 2 AT • (21) 



Note that each one of the four large brackets in(^T|) vanishes separately in the 
(overdetermined) self-dual limit A = (1,1,1,1). In that case, the vanishing 
of each large bracket states the corresponding Bogomol'nyi equation, namely 
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2 tf-h 2 )h r --(l-f 2 )fh = 



1 



[(i-/ a )/»] r + 2(»r , -/» 2 )A = o 



kv 2 - ^ 2 )/], - ^[(i - /W - ^ 2 ) + V 2 h 2 ] = o. 



(22) 
(23) 
(24) 
(25) 



Thus T 44 (|21~1) gives a quantitative measure of the departure of a solution from 



the trivial self-dual field configuration satisfying (|22|)- fl25|).We have plotted 
the profiles of T 44 for each of the four solutions exhibited in Table 1 and 
Figures 1-3, in Figure 4. 



Equations (^)- (]25|) are the spherically symmetric restrictions of the Bo- 
gomol'nyi equations (|[)-(|6]). Specifically, (^2|) pertains to (^), (^) and (^4|) 
to @ and (|25|) to (||). It can be observerved immediately that when each 
of (|22"|)-(|2"5"|) is squared, the sum of the square terms yields the spherically 
symmetric restriction of (|]), namely the one-dimensional system (pi]). The 
corresponding sum a of the cross terms is then the spherically symmetric 
restriction of the topological charge density g = d^Q^ (|3|), divided by 8 ^ 3 . 
Since g is a total divergence, the corresponding one dimensional topological 
charge density which is this quantity a, turns out to be a total derivative 
given by 



° = ~\t ( 3/(1 ~ l f2) - 6(1 " f2)fg2 + (1 " 3/2)/ ^ 4 ) ' (26) 

which is expressed in terms of the rescaled radius p = rjr and the rescaled 
dimensionless function g(p) = rj~ l h(p). The integral / adp, the topological 
charge of the one-dimensional subsystem ( |i0|) , is immediately evaluated using 
the limits (0) and (0) to yield N = 1. 

Having acheived a detailed quantitative understanding of our spherically 
symmetric solutions, we procede to consider the question of inter-instanton 
interactions. 
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2.3 Inter- instant on interaction 



As we saw above, even for the coupling constant vector A = (1, 1, 1, 1), the 
instanton solution is not self-dual and the stress tensor does not vanish. It is 
therefore expected that the force between such instantons is non-vanishing, 
and hence that a dilute gas of such instantons can be constructed. 

In this paper, we shall adapt Polyakov's0 construction of a dilute gas 
of instantons on R3 to our problem on R4. As in Ref.||, we will set out to 
construct a Coulomb gas of intantons. Given that our solution is expressed 
by the Ansatz (||) and obeys boundary conditions (p~3|) and ([[!), we opt to 



work in the Dirac-string gauge introduced in Refs.|5|][[T4]]. In this gauge, the 
Higgs field in the r >> 1 region $(00) = 77757^3^ is gauged to be a constant 
valued field w $ = 777574, by rotating it in the direction of the x 4 axis under 
the action of a suitable 50 (4) gauge rotation uj. As a result the Higgs 
potential and the covariant derivative of the Higgs field vanish, reducing the 
(asymptotic r >> 1) action density ([[]) to the only remaining term Fj^ . 

The gauge field connection and curvature in the Dirac-string gauge are 
calculated from @ with f(r) = according to (|T3|). Just as the 577(2) sym- 
metry of the connection in the Georgi-Glashow model on IR 3 breaks to U(l) 
for the asymptotic fields in the Dirac string gauge, so the SO (A) symmetry 
of the connection @ breaks to 50(3) here. This 50(3) valued asymptotic 
connection and its curvature are given [H]|rj|, respectively, by 



" Ai = (TT^ 7 '^ ' " At = °' (27) 

J F i:j =~ (jij + _L-;E [i7i]A .^ , "F iA = ^JijXj, (28) 



where the same notations as above are used, and with the index /j = i, 4, i — 
1, 2, 3. The Dirac-string singularity along the negative 24 axis is manifest in 
(P7[). Indeed the SO(d) gauge field on for arbitrary d breaks down to an 
SO(d— 1) asymptotic field in the Dirac- string gauge, given by (|28|) in which 
the subscripts 4 are replaced by d, the gamma matrices are understood to be 
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the (/-dimensional gamma matrices and the singularity in the connection is 



on the negative x^-axis. This is explained in detail in Ref.[14 . Note that the 



apparent lack of rotational invariance in (28) is just a gauge artifact, and it 



is easy to check that gauge invariant quantities such as Tr F 2 V , Tr*F fll> F^, 
Tr F {lv F V pF p ^ etc., are SO (A) scalars. Only in the d = 3 case does fl28|) take 
a manifestly SO (3) invariant form, namely 

The components of the curvature in (pB|) are given in Cartesian coordinates 
and are valid in the region r >> 1. It is therefore desirable to express the 
components of the curvature on the 3-sphere S 3 at infinity in polar coordi- 
nates like the Wu-Yang monopole|]l5| whose only non- vanishing component 
on the 2-sphere S 2 at infinity is ^Fg^ = |cr 3 sin 9. In terms of the two polar 
coordinates ifi and 9, and the azimuthal coordinate <fi on S 3 at infinity, the 
non-vanishing components of the curvature ^F are 



"F^g = sin ip(j 31 cos <fi + 723 sin <fi) , (29) 
= - sin 2 ip sin 9 (712 cos 9 + (731 sin + 723 cos <fi) sin 9) , (30) 
u Fij = — sin ip sin 9 (—712 sin 9 + (731 sin phi + 723 cos (f>) cos 9) , (31) 

Since our instantons are exponentially localised, c.f. ( |T7D and (|I8D, we 
shall consider that a gas of non overlapping instantons at positions {x a }, 
with x a b = \ x a ~ x b \ >> V' 1 can be described by the linear superposition 

= ]Tg a u F^(x - x a ), {q a = ±1). (32) 

a 

Following Ref.[0, we argue that the first contribution to the action comes 
from the sum over a of the action integrals of the one-instanton solution. The 
first, dominant contributions to these integrals comes from the integration of 
the action density 

£W = J2Tr{F^ [u (x - x a ), F pa] (x - x a )} 2 , (33) 
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where the integration is restricted to inside the 4-dimensional spheres of radii 
R around each instanton, such that rf x << R << x a &. The second contribu- 
tion comes from the inter-instanton interaction terms in the action density 
calculated from fl32|) , integrated over the 4-dimensional volume outside these 
spheres with radii R. The second contribution to the action is the integral 
of the following density 

Cf ) = Tr{F p[u (x - x a ),F pa] (x - x b )}{F p[u (x - x c ),F pa] (x - x d )}, (34) 

a,b,c,d 

where at most two of the four summation indices a, b, c and d must be dif- 
ferent. Otherwise, when a = b = c = d, the sum is over one index only and 
( |33|) and ([34]) coincide, and the integral in the region outside R is negligible 
compared with that in the region inside R. 



It turns out that when the asymptotic field strengths (28) in the Dirac- 
string gauge are subsituted into ( |34"D the later vanishes exactly, leading to 
vanishing contribution to the action due to inter-instanton interactions. This 
is because the 4-form curvature F pvpa constructed from two distictly situated 
2-form curvatures F pu (x— x a ) and F pa (x— x b ) itself vanishes exactly. This last 
statement can be understood more succinctly by noting that the asymptotic 
curvature 2-form given by ([2~9|)-(|3~1~|) is defined with respect to the three coor- 
dinates {if), 8, 0} on S 3 and to construct a non- vanishing curvature 4-form 
w Fabcd we need at least four coordinates. 



3 The extended models 



The status of the unit charge spherically symmetric instanton solution of 
the basic SO{4) x U{1) model (p]) is that, while the individual lumps are 
exponentially localised permitting the construction of a dilute instanton gas, 
the asymptotic interactions of these instantons are not sufficiently strong to 
lead to a Coulomb gas, in contrast to the situation in the two well known 
models, namely the 0(2) model in two|[J] and the Georgi-Glashow model in 
t/ireep] dimensions. Since it is our eventual aim to construct a Coulomb 
gas of instantons in four dimensions, the basic model ([]]) must be modified. 
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The search for such possible modifications is the purpose of this Section. 
The asymptotic solutions of these will be given in Subsection 3.1 and their 
numerical solutions in Subsection 3.2. 

The modification of (|I[) will take the form of an extension of (jl]). The 
extended model will consist of (|I|) plus some other gauge invariant terms 
depending on (A^, $), such that at least one term in it depends only on the 
curvature and is independent of the Higgs field This is because in 
the Dirac-string gauge (^7|) all gauge invariant quantities depending on $ 
vanish. If the extension consists of a positive definite quantity, then it would 
be expected that the topological lower bound given by the surface integral of 
(H) remains valid and hence that the new instanton will be classically stable. 

To orient ourselves we reconsider the situation in the basic model. There 
the curvature decays with an inverse square power like the monopole|3| on 
R3 and like the latter the asymptotic connection behaves as 

= 2 9 ~ ld » 9 ' 9 = V 7 *" 9 ~ l = W» 

with cr M = (iff, 1), ando^ = (—iff, 1). After passing to the Dirac-string gauge 
the only term in ([!]) which contributes, namely w F tlupcr , decays with the eighth 
power of r and hence the 4-dimensional integral of it decays with the fourth 
power of r. It is clear why this decay is too strong to result in a Coulmb 
potential since the latter is characterised by its inverse square behaviour in 4- 
dimensions. This leads us unambiguously to the criterion that the extension 
to our basic model must include a term decaying with the sixth power of r, 
and since we expect the asymptotic properties of the extended model to be the 
same as those of the basic model (0), i.e. inverse square decaying curvature, 
the extension must include the third power of the curvature. Inevitably this 
means that a new dimensional constant different from the Higgs vacuum 
expectation value r) must be introduced. We shall refer to this as criterion 
(A) below. 

The other criterion, (B), which we shall require is that the density in ques- 
tion involve no higher powers of the velocity fields than the second. In other 
words, for fixed /x and v, no higher powers of (d p A u ) 2 should occur in the new 
density. This criterion is respected by all Skyrme-like models, including the 
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Skyrme model [|I(J as well as the hierarchy of Yang-Mills models and their 
descendents]7]]. In the context of some more phenomenologically motivated 
considerations, this criterion could be relaxed, but this is not the case in the 
present work. The reason for requiring this criterion is that in its absence the 
definition of the canonical momentum-fields would lead to problems in, say 
the collective coordinate quantisation and more specifically in the analysis 
of the fluctuation spectrum when this is eventually performed. While this 
analysis is beyond the scope of the present work, we envisage that it should 
in principle be accessible and hence insist here on this criterion (B). 

Having stated our main criteria (A) and (B) for the extension, we in- 
troduce the other criteria. Obviously this density must be both Lorentz 
invariant, and gauge invariant. In addition it would be an advantage from 
the viewpoint of esablishing the existence of the instanton of the extended 
model, if this density was positive definite by construction. As we shall see 
below, it will not be possible to satisfy this last criterion. 

Subject to the constraints of Lorentz and gauge invariances, criterion(A) 
narrows the choices for candidates down to the following three densities: 



Tr F uX F x ^ v , Tr F^ x F Xu *F^, Tr (D x F flu )(D x F l 



The last of these is equivalent to the non-local effective action density for 
low momentum gluons POfl |21| r TiF llu D 2 F llul while the other possible term 



Tv{D^F^ u ) 2 with the required dimensions is related to the two of the above 
listed terms through 

Trp^) 2 = 2TrF^F uX F Xfl + l -Ti{D x F^) 2 . (35) 

Invoking the criterion (B) eliminates all but the first of these three candidates, 
namely the extension to our basic model (p must be uniquely 



£ 3 = 4TrF^ u F vX F Xfl , (36) 

where K\ is the new constant with the same dimensions as 77. This 
particular dynamical model (B6[) defined on R 6 , has appeared previously 



16 



in the literature [IB] JlTj in a different context where self- dual solutions are 
found [|I6|1||17| and their stability is examined [18]. In particular, the fact that 



criterion (A) is satisfied by Q56J ) was discussed in some detail in Ref . [TBj . In 
the present context on IR 4 , it was first proposed in Ref. 
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Anticipating a result of the next Subsection, namely that the extension 
of (0) by the density (|36f) will result in a power decay for the gauge connec- 
tion function f(r) which in its absence had decayed exponentially according 
to (|T7D, we introduce a second extension to be added to the first one ( j36|) 
with the aim of rectifying this situation. This second extension, which will 
introduce yet another dimensional constant, is 



C 2 = **Tr {(D^Y + X 5 ( V 2 + , (37) 

where k 2 is this second constant with the dimensions of n, and A5 is an 
unimportant dimensionless constant which will be set equal to 1 below. The 



density (p^) is immediately recognised as the usual quadratic Higgs kinetic 
term plus the usual symmetry breaking potential. This term, being entirely 
$ dependent, will have no effect on the inter-instanton interactions and is 
introduced merely to restore exponential decay. Clearly it is not possible 
to add the usual YM term term to (|3~7j) as this will render the action 
logarithmically divergent and would invalidate the instanton solution. 

We note here that the optional criterion of positive definiteness is not 
satisfied by (j36[). This will necessitate some detailed asymptotic and numer- 
ical analysis of the extended systems carried out below, because the lack of 
positive definiteness of the extension means that the topological stability of 
the instanton of the basic model ([I]) does not automatically guarantee the 
same for the extended model. Indeed, if the coupling constant K\ in (j36|) 
is allowed to be large, its lack of posititvity will destroy the stability of the 
instanton in the extended model. Thus it is imperative that this constant k,\ 
be taken to be small. What small means in this context will be determined 
according to quantitative criteria to be applied in Subsection 3.2 where the 
numerical integrations are described, and the extent to which these criteria 
are met will be discussed in detail in Section 4. Here, in justification the 
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extension (|36|) , we suffice to note that the value of the this density remains 
positive as long as the asymtotic conditions fll3l) are respected as we shall 
find to be the case in the Subsection 3.1 immediately below. 

Summarising, we propose two related extensions of ([!]). The first is ob- 
tained by the addition of (|36| ) to ([I]), anmely £/ = £4 + £3, and the second 
by the addition of ( p6|) and (|37f) to ([J), namely £// = £4 + £3 + £2. 



3.1 Asymptotic solutions 

The one-dimensional subsystems L 3 and L 2 arising from the imposition of 
spherical symmetry on the extensions £3 and £2, analogous to the one- 
dimensional subsystem ( p~0| ) of ([!]) are, respectively 



L 3 = ^(l-f 2 )(3f^ + -Jl-fr) (38) 



i 2 = i^r 3 hi + ^fh 1 + 2A 5 („ 2 -fc 2 ) 2 . (39) 



It is important to note that the one dimensional density L 3 remains posi- 
tive provided that the asymptotic conditions (p~3|) are satisfied by the solution, 
inspite of not being a positive definite quantity by construction. We shall 
find this to be the case below when we solve the Euler-Lagrange equations 
5fLi = and SjLji = in the region r >> 1. This will be the justification 
of employing (|36|) as an extension. 



The Euler-Lagrange equations 5fL 3 = for (|3"8"D, and, <5/L 2 = and 
ShL 2 = for (B^), are given respectively by 



S f L 3 = ((1 _ f)f rr _ f f r _ i(i _ /2 )/r + J_(i _ (40) 

5 f L 2 = 3«£r /i 2 / (41) 
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5 h L 2 = n 4 r fefh - 4A 5 r 2 (r? 2 - h 2 ) h - 3rh r - r 2 h„) . (42) 



To start with we dispose of the question of the asymptotic solutions in the 
r << 1 region, to both the extended models Ci and Cn. We have verified 
that the additional terms do not change the behaviours (|1^) and ( pD| ) for 
the functions f(r) and h(r) respectively. There then just remains to find the 
asymptotic behaviours in the r >> 1 region. 



Consider first the extended system Ci = £ 4 + £ 3 . The 5hLj = ^L 4 = 
equations of its one-dimensional subsystem are unchanged and are given by 
fljjD, while the SfLi = SfL^ + 8fL 3 = equations are now given by (pTT]) and 
Linearising the latter around the asymptotic value of f(r) = + F(r) 



we find 

r 2 Frr _ rFr+ ^l_(W^ F = 0j (43) 

which has a power decay solution yielding 

f{r) ~ r 1 -^^ (44) 



provided that a/A7^ > 1- This is not physically unreasonable since the 
constant r\ is expected to be large, while K\ can be taken to be small. In 
practice we have set r] = 1 in all our numerical computations and, taken Ai 
to be of the order of unity, while K\ must be taken to be sustantially smaller 
than unity if the lack of positive definiteness of £3 (|36|) is not to predjudice 
the instanton solution of the basic model. Thus, for the parameters employed 
in our numerical computations, the conditions for (|44]) to imply a rapid decay 
are met. 



Notwithstanding the fact that we expect a rapid decay (0) for the func- 
tion f(r), it would be desirable to extend the model in such a way that like 
the function h(r), the function f(r) also decays exponentially. This is be- 
cause the validity of the approximations that will be made in an eventual 
application of the instanton solutions found here to the construction of a 
dilute instanton gas, rely on the strong localisation of the instantons and the 
best way of meeting this requirement is by exponential localisation. 
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To this end, we consider the extension Cu = £4 + £3 + £2. In this 
case, the one-dimensional Euler-Lagrange equations SfLu = SfL^ + 5fL 3 + 
5 f L 2 = and 8 h L n = 5 h L 4 + 5 h L 2 = are given by ([□]), ([Lj), (|0]) and 
fl4"2"|). Linearising these around their asymptotic values f(r) = + F(r) and 
h(r) = 77 + H(r) we find 

p 2 F pp - pF p -p 4 F = (45) 
<j 2 H aa + 3a# CT - a 2 i/ = (46) 

where we have used the dimensionless rescalings p = ( ) 1 ? ? r anc ^ °" = 

2^j2X^rjr of the radial variable r. Note that the notations p and cr in (^5|) 
and (|46| ) are different from those in (|15|) and (|16|). From the fourth power of 
p in the last term of (|4"5| ) it is clear that this equation cannot be brought to 
the form of a modified Bessel equation so we evaluate its asymptotic solution 
directly, yielding f(r) in the r >> 1 region to be 

K 2 

2 y,2 



f{r) = e 2 ^i" ' , (47) 

which is clearly exponentially localised. Equation (^) on the other hand 
can be brought to the form of a Bessel equation with the solution H(a) = 
(j^Ki^a) yielding 



Hr) -^ 1 -2^r K ^ 2Xmr) )- <48) 

In summary, we see that in the extended model Ci = £4 +£3, the function 
h(r) is exponentially localised according to ( |i~8|) , while the function f(r) is 
power localised according to (|44"|). In the extended model £77 = £4 + £ 3 + £2 
both functions f(r) and h(r) are exponentially localised according to ([47|) 
and (^|) respectively. 

We see that the function h(r) is exponentially localised according to (|i~8|) 



and (^8f) for both the extended models £/ = £4 + £ 3 and Cu = £4 + £3 + £2, 



respectively, while the fuction f(r) is power localised according to (f44|) for 
the model £/ = £4 + £3 and exponentially localised according to ([4~7|) for 
£// = £4 + £3 + £2- 
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Perhaps the most important result of this Subsection is, that the asymp- 
totic behaviours of the function /(r) in both extended models given respec- 
tively by flI§D,(fPD, and, (|rUD,(pETD, ensure that the one dimensional density 
(|38D does not become negative. 



3.2 Numerical solutions 

For the purposes of numerical integrations we have employed the following 
two values of the coupling constant Ki, k1 = 0.1 and 0.01 in(^), which result 
in appropriately small quantitative extensions of the basic model. 

For the extended model £j = £4 + £3 given by ([1]) and fl36|), we have 
integrated equations 5fL I = SfL^ + SfL 3 = and b^Lj = b^L^ = given re- 
spectively by flnp, P0|), and by flT2]), numerically, starting with the functions 
( |19"D and (|20|) near r = and choosing the numerical values of the constants 
A and B so that f(r) and h(r) reach their asymptotic values given by the 
second members of (|T3| ) and (0). We have performed this integration only 
for the dimensionless coupling constants = (1,1,1,1). The numerical 
values of the pair of constants {^4, B} are fixed by the numerical integra- 
tions, and are listed in the first and second colums of Table 2. The profiles of 
the functions f(r) for this extended model are not given since they are quali- 
tatively the same as those in Figure 1, while the function h(r) is the same as 
that for the basic model and is given by the A = (1, 1, 1, 1) curve in Figure 2. 
The numerical values of the total actions Si(ki = 0.1) = and Sj(ki = 0.1) = 
pertaining to these two solutions are listed in the third column of Table 2. 

For the extended model Cu = £4 + £3+ £2 given by (JTJ), ( |5BD and (|5TD, we 
have employed the same two values of the coupling constant K\ — 0.1, 0.01 as 
above and for each of these the values of the coupling constant k 2 = 0.25 and 
k 2 = 0.5. We have integrated the two second order equations SfLjj = 5/L 4 + 
6fL 3 + 5 f L 2 = given by (0) (0) and ©, and, 5 h L n = 6 h L 4 + 5 2 L 2 = 
given by ( |I2"D and (^). Again we have performed this integration only for 
A = (1,1,1,1), and the values of the pair of numerical constants {A, B} 
are fixed by the numerical integration and are listed in the first and second 
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columns of Table 3. Again the profiles of the functions f(r) and h(r) are not 
exhibited since they are qualitatively the same as those in Figures 1. and 2. 
The numerical values of the total actions pertaining to these two solutions 
are Sh{k\ = 0.1) = and Sh(ki = 0.01) =. 

Since the density £3 ( pop is not positive definite, the stability of the 
instantons of both the models decsribed by the densities £1 and £n are 
not guaranteed by the stability of the instanton of the basic model. We 
argued above however, that if the coupling constant K\ is small enough, the 
instantons of the extended models will nevertheless be stable. To this end, 
the values of the total actions corresponding to each instanton given above 
are relevant. Quantitatively, we seek to demonstrate that the addition of 
the extension densities (^) and ( |5?D to the basic density ([[[) does not result 
in an appreciable change in the value of the action integrals. The detailed 
comparisons are deferred to Section 4. 

Before preceding to Section 4 however, we must perform one further 
numerical operation. From the viewpoint of these quantitative comparisons, 
the difference of the extended model £j from the basic model is just the non- 
positive density £3 and hence the interesting quantity is the difference of their 
action integrals. The extended model £n however is the result of the addition 
of the non-positive density £3 to the positive definite density £m = £4 + £2 
and not to the basic density alone. The interesting quantity in this case 
therefore is the difference of the action integral of £u from that of the system 
described by £ni, which we have not studied numerically hitherto. This we 
now do by integrating the Euler-Lagrange equations SfL in = 5fL i +5fL 2 = 
given by ( |TTD and fl4"T|), and ShLjn = b^L^ + 5hL 2 = given by ( |Ti2"D and 
(f42D . We have performed the numerical integrations for the value of = 
(1, 1, 1, 1) only, and for the values of the coupling constant n 2 = 0.25 and 
k 2 = 0.5. The numerical values of the constants {^4, B} for these solutions are 
listed in the first and secon columns of Table 4, while the values of the total 
action arelisted in the third column of Table 4. We do not plot the profiles 
of the functions /(r) and h(r) in this case because they are qualitatively the 
same as those in Figures 1 and 2. 
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4 Summary and discussion 



We have made a detailed numerical study of unit topological charge instanton 
solutions of the 5*0(4) x U(l) Higgs model|| and its extensions |T9[1 on IR 4 . 
The extensions to the basic model are devised such that the dynamics of the 
extended model would in principle be capable of supporting a dilute Coulomb 
gas of interacting instantons. This construction however is deferred to future 
work, and in the present work we have restricted ourselves to establishing 
the existence of the classical instanton solutions using numerical methods. 



In Section 2, we have studied the spherically symmetric instanton solution 
of the basic models (|l|) parametrised by a set of dimensionless constants A = 
(Ai, A 2 , A3, A 4 ), for four different such sets \,i = 1,2,3,4 listed in Table 1. 
The results are exhibited in Table 1 and Figures 1 and 2. The unusual feature 
of the system (M is that for the case i = 3, namely for A3 = (1, 1, 1, 1), where 
the lower bound on the action is saturated by the Bogomol'nyi equations (||)- 
(||), the only solution is the trivial constant asymptotic values of the fields. 
That equations are overdetermined can clearly be seen by inspection of 

their spherically symmetric restrictions (|22|) -(25), which are four constraints 
on two functions. This feature is shared with the usual Skyrme model|Tn], in 
which case however the departure of the soliton solution from the self-dual 
field configuration is fixed by the absolute scale there. Here by contrast this 
departure is depends on the dimensionless coupling constant A,, and in this 
sense our model is more akin to the Abelian Higgs model on R 2 which even 
though it does support self- dual solutions, the departure from self-duality 
depends on the value of the dimensionless coupling strength A multiplying 
the Higgs potential. A in that case is the analogue of Aj here. In both models 
the value of the action integral increases with an increase in the value of these 
dimensionless coupling strength(s). This is seen to be the case from the third 
column of Table 1, where all the action integrals S(Xi) are larger than the 
topological lower bounds. In particular the action integral S^A/i) is larger 
than the lower bound 1 in that case, while the actions S(Xi) for i — 1,2,3 
are each larger than their respective lower bounds (Aj) min . To demonstrate 
more quantitatively the departute from the self dual configuration, we have 
plotted the profiles of the one-dimensional restriction of the stress-tensor T44 
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for each solution in Figure 4. 

In Section 3 we have given some arguments in favour of our choice of 
extended models and preceded to study these numerically. It turns out that 
our choice for an extension of ([!]) is narrowed down to the cubic density 
(j3"6|) . This density being non positive definite, there is no guarantee that 
its addition to (|TJ) will not invalidate the topological lower bound. In that 
case the extended model would not support an instanton solution. We have 
argued that for small enough coupling of this term, the instanton solution 
persists in the extended model is the dynamics of the basic model is robust 
enough. The numerical data pertaining to this extended model Ci is given 
in Table 2. We can see that the action of £/ does not differ appreciably from 
the action of the basic model, by comparing the actions Si(X^, K\) in Table 
2 with the action S(\s) in Table 1. Indeed, this difference is smaller for the 
case with the smaller value of the coupling constant k,\, thus justifying our 
claim that for small enough modifications of the basic model we can expect 
to have intanton solutions. 

The function h(r) parametrising the instanton of the extended model £/ 
is localised exponentially according to (|T8|), but the function /(r) is power 
localised according to (fBl) . This has motivated us to consider yet another 
extension, for which both functions h(r) and f(r) are both exponentially 
localised. This is the model Cu which results from the further addition of 
the density (|37]) to Ci and whose exponentially decaying solutions are given 
by (f47|) and ([Kf). The numerical data for this model is given in Table 3. To 
carry out the corresponding quantitative checks as in the previous case, we 
need to compare the values of the actions Sjj(X 3 , k,i) with the action integral 
of the positive definite density defined by the sum of ([D and the positive 
definite density £2 defined by fl3"T|). We have denoted this system by Cm 
and listed the numerical data pertaining to it in Table 4. We must therefore 
compare the actions Su(\ 3 ,k,i) in Table 3 with the corresponding actions 
Sin{\3, K2) in Table 4. The qualitative conclusions are exactly the same as 
those for the extended model Cj, namely that the action in the extended 
model Ln is close enough to that of the basic model, that £77 also supports 
instanton solutions. 



24 



Quantitatively, the differences Sj ( A3, = 0.01, k 2 = 0) — S(\s,Ki = 
0,k 2 = 0) = 67 x 10~ 6 and S n {X 3 ,Ki = 0.01, k 2 = 0.25) - S///(A 3 ,/ci = 
0, k 2 = 0.25) = 68 x 10~ 6 are practically equal, implying that the system 
([!]) is equally robust against the addition of (|36D as against the addition of 



( |36"D and fl3"T|), even though in the second case the exponential decay was 
better satisfied. Evidently, this means that the value chosen for the coupling 
constant Hi was sufficiently small so that the extended systems arising from 
the addition of the non-positive definite F 3 term (^) result in new systems 
that also support instanton solutions. To test the validity of this criterion, 
consider the differences of the corresponding actions for the extended models 
with a larger value of the coupling constant K\. Thus 57 (A3, /«i = 0.1, k 2 = 
0) - S(\ 3 , ki = 0, k 2 = 0) = 6665 x 10~ 6 and S//(A 3 , k x = 0.1, k 2 = 0.25) - 
Sju(X 3 ,K,i = 0, k 2 = 0.25) = 6685 x 10~ 6 , which are again nearly equal to 
each other but are very considerably larger than the former pair of differences. 
This completes our justification of the extended models Ci = £4 + £3 and 
Cn = £4 + £3 + £ 2 . 

Having established that the extended models also can support instantons 
with the same asymptotic properties as the instantons of the basic model, 
the remaining task is to compute the volume integral 



Sint = J d A x Tr ^( x ~ x a)F v \{x - x - x b )F x ^{x - x c ) (49) 



in the Dirac string gauge. Clearly Si n t in fl49|) is the only nonvanishing 
contribution to the action due to inter-instanton interactions since in the 
Dirac string gauge the asymptotic contributions of all Higgs dependent terms 
vanish and the density of the only other Higgs independent term in either 
extended action density coming from F^ vpa also vanishes according to our 
arguments at the end of Subsection 3.2, namely (0). The integrand in 
(H|), which is readily calculated by substituting (|28|) in ([49]), is a quantity of 



considerable complexity and since we will not perfom the integration here we 
do not write it down. If one set any two of the summation indices a, b, c, in 
( f4"9"D the same, namely restricted to "two neighbour" interactions, then the 
integrand would become singular at the origin in the usual way. Thus the 
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interaction term will have to take account of "three neighbour" interactions. 
This situation may be changed by replacing the cubic extension term (|36| ) in 
(pE9D by the alternative density F^ U D 2 F^ U appearing in Q55D 

Sint = J d A x ^D x F llu (x - x a )D x F llu (x -x-x b ), (50) 

a, b 

which is a "two neighbour" interaction. The details of these contributions 
to the action as well as their consequences in the resulting effective theory of 
Strong interactions will be investigated and given elsewhere. 

AcknowledgementsWe thank J. Burzlaff for his cooperation at the early 
stages of this work. 



26 



References 

[1] A.A. Belavin, A.M. Polyakov, A.S. Schwarz and Yu.S. Tyupkin, Phys. 
Lett. B59 (1975) 85. 

[2] A.M. Polyakov, Nucl. Phys. 120 (1977) 249. 

[3] G. 'tHooft, Nucl. Phys.79 (1974) 276; A.M. Polyakov, JETP Lett. 20 
(1974) 194. 

[4] V.L. Berezinsky, JETP 30 (1970) 493; J.M. Kosterlitz and D.J. Thou- 
less, J. Phys. C6 (1973) L97. 

[5] G.M. O'Brien and D.H. Tchrakian, Mod. Phys. Lett. A4 (1989) 1389. 

[6] D.H. Tchrakian, "Yang-Mills Hierarchy", Int. J. Mod. Phys. A (Proc. 
Suppl.) 3A (1993) 584, and references therein. 

[7] D.H. Tchrakian, " Skyrme-like Models in Gauge Theory", in "Constraint 
Theory and Quantisation Methods", eds. F. Colomo, L. Lusanna and 
G. Marmo, World Scientific, Singapore, 1994, and references therein. 

[8] T.N. Sherry and D.H. Tchrakian, Phys. Lett. B 295 (1992) 237. 

[9] L. Jacobs and C. Rebbi, Phys. Rev. B19 (1979) 4486. 

[10] T.H.R. Skyrme, Nucl. Phys. 31 (1962) 556. 

[11] D.H. Tchrakian and A. Chakrabarti, J. Math. Phys. 32 (1991) 2532. 

[12] R. Jackiw and C. Rebbi, Phys. Rev. Lett. 37 (1976) 172; C.G. Callan, 
R.F. Dashen and D.J. Gross, Phys. Lett. B63 334. 

[13] see A.M. Polyakov, "Gauge Fields and Strings", Harwood Academic 
Publishers, Chur, 1987. 

[14] Zh-Q. Ma and D.H. Tchrakian, Lett. math. Phys. 26 (1992) 111. 

[15] T.T. Wu and C.N. Yang, J. Math. Phys. 15 (1974) 53. 

[16] C. Saclioglu, Nucl. Phys. B227 (1986) 487. 



27 



[17] K. Fujii, Lett. math. Phys. 12 (1986) 363. 

[18] Zh-Q. Ma and D.H. Tchrakian, Lett. Math. Phys. 19 237. 

[19] G.M. O'Brien and D.H. Tchrakian, "Exponentially localised instantons 
in a hierarchy of Higgs models" , in " Geometry of Constrained Dynamical 
Systems", ed. J.M. Charap, Cambridge University Press, Cambridge, 
1995. 

[20] D.I. Diakonov and V. Yu. Petrov, Phys. Lett. B242 (1990) 425. 

[21] M. Baker, J.S. Ball and F. Zachariasen, Nucl. Phys. B229 (1983) 445. 



28 



TABLES: 



Table 1 



K 


A 


B 


S(Xi) 


0.5 


0.300044 


0.653908 


1.023269 


0.8 


0.324460 


0.650557 


1.611897 


1.0 


0.335691 


0.648918 


2.002486 


1.2 


0.344599 


0.647561 


2.392250 



Table 2 



K 2 =0 


«1 




B 


Sj(A 3 ,ki) 


0.01 


0.335669 


0.648919 


2.002553 


0.1 


0.333488 


0.648989 


2.009151 


Table 3 


k 2 = 0.25 




A 


B 


Sn(X 3 , «l) 


0.01 


0.336761 


0.651287 


2.010352 


0.1 


0.334577 


0.651353 


2.016969 


«i = 0.50 




A 


B 


<S//(A 3 , Kl) 


0.01 


0.350511 


0.679313 


2.108078 


0.1 


0.348347 


0.679341 


2.115000 



Table 4 



Kl = 








Slll(\3, K2) 


0.25 
0.50 


0.336783 
0.350533 


0.651286 
0.679313 


2.010284 
2.108006 



29 



FIGURE CAPTIONS: 



Figure 1. Profiles of the function f(r) for the systems characterised by Ai = 
(0.5, 0.5, 0.5, 0.5), A 2 = (0.8, 0.8, 0.8, 0.8), A 3 = (1, 1, 1, 1), A 3 = (1.2, 1.2, 1.2, 1.2) 
respectively from right to left. 

Figure 2. Profiles of the function h(r) for the systems characterised by Ai = 
(0.5, 0.5, 0.5, 0.5), A 2 = (0.8, 0.8, 0.8, 0.8), A 3 = (1, 1, 1, 1), A 3 = (1.2, 1.2, 1.2, 1.2). 
All curves strongly overlapping. 

Figure 3. Profiles of the action densities of the solutions to the systems char- 
acterised by Aj = (0.5,0.5,0.5,0.5), A 2 = (0.8,0.8,0.8,0.8), A 3 = (1,1,1,1), 
A 3 = (1.2, 1.2, 1.2, 1.2) respectively in order of the heights of the curves. 

Figure 4. Profiles of T44 given by fl2~l|) for the systems characterised by Ai = 
(0.5, 0.5, 0.5, 0.5), A 2 = (0.8, 0.8, 0.8, 0.8), A 3 = (1, 1, 1, 1), A 3 = (1.2, 1.2, 1.2, 1.2) 
respectively in order of the lowest to highest curves. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-th/9603047vl 



